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Abstract 

Closed formulas are derived for the field in the focal region of a diffraction limited lens, such that 
the electric field component in a given direction at the focal point is larger than that of all other 
focused fields with the same power in the entrance pupil of the lens. Furthermore, closed formulas 
are derived for the corresponding optimum field distribution in the lens pupil. Focused fields with 
maximum longitudinal or maximum transverse are considered in detail. The latter field is similar, 
but not identical, to the focused linearly polarized plane wave. 

PACS numbers: 42.79.Jb, 41.85. Ct, 42.25Hz, 42.25.Ja, 42.79.Ci 
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I. INTRODUCTION 



When a linearly polarized plane wave is focused by a diffraction-limited lens, the intensity 
distribution in the focal plane is in the scalar theory the well-known Airy pattern. However, 
when the lens has high numerical aperture, the rotation of polarization must be accounted 
for and the vector diffraction theory of Ignatowsky and Richards and Wolf has 

to be applied to obtain the field distribution in the focal region. We then get three electric 
and three magnetic field components in the focal region. When the beam in the lens aperture 
is uniformly linearly polarized plane wave, the dominant electric field component in the focal 
region is found to be parallel to the polarization direction of the incident plane wave. But, 
as the numerical aperture increases, the maximum value of the longitudinal component of 
the electric field in the focal plane becomes quite substantial, although it vanishes at the 
focal point itself. 

An appropriately shaped focused spot is essential in many applications such as optical 
recording, photolithography and microscopy. Furthermore, a field in focus with maximum 
electric component in a specific direction is important for manipulating single molecules 
and particles, and in materials processing j^, p, Ql, js], 9]. The focused wavefront can 
be tailored by setting a proper amplitude, phase and polarization distributions in the pupil 
of the focusing lens. Nowadays it is possible to realize almost any complex transmission 
unction in the pupil plane, using for example liquid crystal-based devices [10|, [ll|, |12| . 

The optimization of focused fields has been studied by other authors. 

In this paper, we maximize a specified component of the electric field in the focal point 
of a diffraction limited lens. First we consider fields in free space or in homogeneous matter, 
without taking into account the way these fields are realized, in particular without consider- 
ing the lens. We merely suppose that, with respect to a cartesian coordinate system (x, y, z), 
the fields considered consist of plane waves propagating in the positive z-direction and have 
wave vectors with angles with the positive 2;-axis that do not exceed a specific maximum 
angle, i.e., the numerical aperture of the plane waves is restricted. The optimization prob- 
lem is then to find the complex plane wave amplitudes such that, for a given direction in 
space and for given mean flow of electromagnetic power through a plane 2;=constant, the 
amplitude at some point (chosen at the origin) of the electric field component that is parallel 



2 



to the chosen direction is larger than that of any other field for the same numerical aperture 
and the same mean total power flow. We shall derive closed formulas for the plane wave 
amplitudes of the optimum fleld. 

The solutions of the optimization problems for the fleld propagating in homogeneous space 
is rigorous results since they are derived from Maxwell's equations without any further as- 
sumptions. Next we will consider the realization of the optimum fleld using a diffraction 
limited lens with the origin as focal point. By using the vector diffraction theory of Igna- 
towksy and Richards and Wolf, closed expressions for the optimum pupil distributions will 
be derived, which after focussing give the optimum fleld component in the focal point of 
the lens. In contrast with the solution in terms of the plane wave expansion, the formulas 
for the optimum pupil flelds are approximate since they are based on the vector diffraction 
theory which is an approximate theory that is valid for lenses of which the focal distance 
and the pupil radius are many wavelengths. 

When one considers the focussing by a lens it is obvious that the plane wave expansion 
in image space has flnite numerical aperture of that of the lens: NA = rasinamax, where 
n is the refractive index in image space and Omax is half the top angle of the cone with 
top the focal point and base the pupil. But also when one would consider only waves in 
free space without a focussing lens, there is a good reason to restrict the flelds to flnite 
numerical aperture, in that case in particular to NA = n. In fact, when NA > n a, part of 
the evanescent waves are taken into account in the expansion. One can then construct flelds 
with a given power which have arbitrary large components. The evanescent waves do not 
contribute to the total power and hence one can increase their amplitudes by any desired 
amount to achieve arbitrarily high local flelds. Stated differently, by constructing suitable 
time-harmonic source distributions that emit singular flelds with flnite power flow, one can 
achieve arbitrary large fleld components by approaching these source distributions. At small 
distances to the sources the evanescent waves play of course a major role. 

Among the directions for the optimized electric fleld component, two are of particular 
interest, namely the directions parallel and perpendicular to the optical 2;-axis. These di- 
rections are also called the longitudinal and transverse directions. The fleld with maximum 
longitudinal component has been discussed in [l^ , but without derivation. In this paper de- 
tails of the derivation are provided and the optimization problem is generalized to arbitrary 
directions of the electric fleld vector. 
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As was announced in 



15| , the pupil field that when focused gives maximum longitudinal 



component, is radially polarized. This means that in all points of the lens pupil, the electric 
field is linearly polarized with the electric field pointing in the radial direction. Furthermore, 
the electric fields in all points of the pupil are in phase and the electric field amplitudes are 
rotationally symmetric. The amplitude of the electric pupil field vanishes at the center of 
the pupil and is a monotonically increasing function of the radial coordinate. The shape of 
this function depends on the numerical aperture. 

It wa. noted by several autho. Q. B. Q. Q that when a .adiaUy polan.ed bean, ,s 
focused, the distribution of the longitudinal component can be considerably narrower than 
the focussed spot obtained by focusing a linearly polarized plane wave. With the develop- 
ment of a new generation of photoresists j20|, it is possible to control the photosensitive 
material in such a way that it will react to only one of the polarization components of the 
electric field. Materials with molecules having fixed absorption dipole moments have been 
applied in to be able to probe field components individually. When this component is 
the longitudinal component, a tighter spot can thus be obtained than with the classical Airy 
pattern. 

Often the amplitude distribution of the radially polarized beam in the pupil plane is 
chosen to be a doughnut shape or a ring mask function 17| . But these distributions do not 
give the maximum possible longitudinal electric field component in focus for given power 
and its amplitude as a function of the radial pupil coordinate differs from the optimum 
function derived in the present paper. 

The other case of particular interest is the optimization of the transverse electric field 
vector. Since the optical system is assumed to be rotationally symmetric around the optical 
axis, we may choose this direction parallel to the x-axis. The solution of the optimization 
problem is then the field for which the amplitude of the x-component of the electric field 
in the focal point is maximum for the given numerical aperture and the given total power 
flow. We will show that the corresponding pupil field is linearly polarized with direction 
of polarization predominantly, although not exactly, parallel to the x-axis. Therefore, the 
focused optimum field is similar to the vectorial Airy pattern of a focused linearly polarized 
plane wave, although it is not identical to it. 

The paper is organized as follows. In Section 2, we will formulate the optimization 
problem and we will prove that the optimization problem has one and only one solution. 
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In Section 3, we will apply the Lagrange multiplier rule to obtain closed formulas for the 
plane wave amplitudes of the optimum field and for the optimum field distributions near 
focus. In Section 4, we will study the optimum fields in the focal region, in particular their 
mean energy flow. Then, in Section 5, we apply the vector diffraction theory of Ignatowsky, 
Richards and Wolf to derive the electric field distribution in the pupil of the lens that, when 
focussed, yields the maximum field component in the focal region. 



II. FORMULATION OF THE OPTIMIZATION PROBLEM FOR ARBITRARY 
ELECTRIC FIELD COMPONENT 

We begin with some notations. Consider a time-harmonic electromagnetic field in a 
homogeneous unbounded medium with real refractive index n (i.e., the material does not 
absorb electromagnetic radiation of the given frequency): 

S{v,t) = Re [E(r)e-*^*] , (1) 
n{r,t) = Re [H(r)e"*^*] , (2) 

where a; > 0. As stated in the Introduction, the lens is first not considered in the op- 
timization problem. It is merely assumed that, with respect to the cartesian coordinate 
system (x, y, z) with unit vectors x, y and z, the electromagnetic field ([T]), ^ has numerical 
aperture NA < n and that the plane wave vectors have positive ^-component: 

^("■^ = i i A{k,,ky)e'''-''dhdky, (3) 

H(r) = / / kxA{k,,ky)e'''''dk,dky, (4) 

47r CU/io J yyfc|Tfef<fconsinQ„,ax 

where k = {kx, ky, kz) with 

kz = {kln^ — fc^ — kyY^"^, (5) 

with ko = Uy/eQfiQ = 2ti/\q where Aq is the wavelength in vacuum, and where NA = 
nsinamax with Omax the maximum angle that the wave vectors make with the positive z- 
direction. If NA = n we have Omax = vr/2 and the plane wave spectrum then consists of 
all homogeneous plane waves that propagate in the non-negative ^-direction (there are no 
evanescent waves in the expansion). When NA < n, the cone of allowed wave vectors has 
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top angle Omax = arcsin(NA/n) < 90°. Because the electric field is free of divergence we 
have that 

A ■ k = 0. (6) 

We shall use spherical coordinates in reciprocal k-space: 

k = sinacos/3x + sinasin/5y + cosaz, (7) 
a = cos a cos/5x + cosa sin/5 y — sin a z, (8) 
j3 = — sin /3 X + cos j3 y, (9) 

where < a < Omax and < /5 < 27r are the polar and azimuthal angles, respectively. 
Conversely, we have 

X = sin a cos /3 k + cos a cos /3 CK — sin /3 (10) 
y = sin a sin /5 k + cos a sin /5 q: + cos /5 /3, (11) 
z = cos a k — sin a a. (12) 

Note that {k, a, (3} is a positively oriented orthonormal basis: 

k X a. = ^, a X ^ = 'k, ^ x 'k = a. (13) 

Furthermore, k = konh and the Jacobian of the transformation {a, (3) ^ {k^, ky) is: 

cos a cos f3 — sin a sin f3 
cos a sin /3 sin a cos /3 

so that 

dkx dky = k^v? cos a sin a da djS. (15) 

Because of (E]) we have: 

A{a,l3) = Ao,{a,l3)a + Ap{a,l3)p, (16) 
for some functions and Aj^. Then, using (fT3l) : 

kxA = kQn{-ApdL + Aaj3). (17) 

The plane wave expansion can thus be written as 

E(r) = T2 / / {AadL + Afi^) COS asm. a e'^''' da d[5, (18) 

H(r) = -2 — / / {-Apa + A^p) cos a sin a e^-"" da d(3. (19) 





dkx \ 




j da 

1 dky 


a/3 ' 

dky 1 




\ da 


a/3 / 





(14) 



The a-component is parallel to the plane through the wave vector and the 2;-axis, whereas 
the /3-component is perpendicular to this plane. 

Let V = f + VyY + f 2Z be a real unit vector. We consider the projection of the electric 
field at the origin at time t = on the direction of v: 

E(0) ■ V = / / [Aa{a, (3) Va + Ai3{a, (3) vp] COS a sin a da d(3, (20) 



^0 "'0 ^0 



where 



Va = V ■ a. = Vx cosa COS P + VyCosasin P — VzSina, (21) 
= V ■ f3 = —Vx sin (3 + Vy cos [3. (22) 

We will consider E(0) ■ v as a (linear) functional of A = Aadt + Apf5, which for brevity we 
will denote by F{A). Hence, 

F{A) = -To \ / [Aa{a -I (3) Va + Ai3{a,(3)vi3\ COS a sin a da d(3. (23) 
-^0 Jo Jo 

Next, we calculate the total mean fiow of power through a plane z =constant. The total 
mean power fiow is obtained by integrating the normal component of the vector (l/2)ReS 
over the plane z = constant, where S = E x H* is the complex Poynting vector. By using 
Plancherel's formula, the integral of Re S over this plane can be written as an integral over 

and Ivy'. 



/ -Re [S(r)] dxdy = -Re / / E(r) x B.{r)* dxdy 

■ooJ —00 2 2 J —ooJ —00 



\Re / / A{kx,ky)e''^' X 



X A(^kx, ky) 



UJflQ 



Q 2 / / 

= l-^ff \Aikx,ky)\^kdkxdky, (24) 

where we used that k is real and A{kx, ky) ■ k = 0. The total time-averaged fiow of energy 
in the positive 2;-direction through the plane z = constant is given by the z-component of 



/ / l-Re[Sz{r)]dxdy = \A{kx, ky)\^ k^dkxdky 

J —00 J —00 



3 / \ /-Qmax /-Ztt 

[|y4Q,(a, /5)p + |A^(a, /3)p] cos^ a sin a da d(3. 



_ ( £!I 
2\l V/io 



JO 



(25) 



This is independent of the plane z = constant, as should be in a medium without losses. 

The quantity F{A) = E(0) ■ v is the complex electric field component in the direction of 
V at time t = 0. Without restricting the generality we may assume that F{A) is real. If it 
were not real, a time shift could be applied to make it real. Hence we may assume that 

Im [E(0) ■ v] = / / Im [Aa{a, (3)va + Ap{a, (3)vi3] cos a sin a da d(3 = 0. (26) 
Jo Jo 

The optimization problem is to find the plane wave amplitudes A = AaCt + Aj^fB for which 
the electric field component at the origin that is parallel to the direction of v is larger 
than for any other field with the same mean power fiow through a plane z = constant and 
the same numerical aperture. To formulate this problem mathematically, we introduce the 
space Ti of plane wave amplitudes A = AaCt + ApfB which have finite mean fiow of power 
through the planes z = constant: 

H = {A = AaO. + A/sP; / / [\Aa{a,/3)\'^ + \A/3{a,(3)\'^] cos"^ a sin a da d(3 < oo}, 

Jo Jo 

(27) 

and we define Tio as the subspace of Ti. consisting of all A which satisfy fl26l) . Then 7^ is a 
Hilbert space with scalar product 

< A, B >7^= / / [Aa{a, [3) Ba{a, (3)* + Af3{a, (3) Bi3{a,/3)*] cos"^ a sinadad(3,{28) 
Jo Jo 

and Tio is a closed subspace of Ti. Note that the constraint (126!) means that Im(A) is 
perpendicular to the vector field v/ cos a in the space Ti, i.e. Im(A) is perpendicular to 
v/ cos a in the sense of scalar product ( !28l) . 
Define the quadratic functional 

r\ f r?^ / f \ /-"max r2TT 

P(A) = — 2 — / / [\A^{a,f3)\^ + \A^{a,/3)\^] cos^asinadadp, (29) 
2Ao \fioJ Jo Jo 

which is the mean power fiowing through a plane z = constant for fields with plane wave 
amplitudes A = AaOc + A/^(3. Then the optimization problem is to find, for given Pq > 0, 
the solution of 

(*) max F(A), under the constraint -P(A) < Pq- 

AgWo 

For any solution of problem (*) the equality -P(A) = Pq holds, because otherwise A could be 
multiplied by the number (Po/-P(A))^/^ > 1 and this would increase the value of F without 
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violating the constraint on the energy. Hence it does not matter whether we impose the 
equahty constraint -P(A) = Pq or the inequahty constraint -P(A) < Pq on the mean flow of 
energy. 

It is not completely obvious that problem (*) has a unique solution since it is posed in 
a linear space Tio of infinite dimension. However, there is a functional analytic theorem 
which states that a continuous real linear functional attains its supremum on a sphere in 



a Hilbert space and that the solution is unique [2]J. Since the functional F is linear, real 
and continuous with respect to the norm on Ti, and since the feasible set of problem (*) is 
a sphere in Tio, this theorem applies to our problem. Hence the optimization problem has a 
unique solution. In the next section we shall compute the solution. 

III. OPTIMUM PLANE WAVE AMPLITUDES 

Since F is a linear functional, the Frechet derivative of F at A in the direction of B is 
simply -F(B), i.e. 

5F(A)(B) = F(B) = -2 / / [Ba{a,/3)vc, + Bp{a,(3)v/3] cosasinadadp. (30) 

■^0 Jo Jo 

The Frechet derivative of the quadratic functional -P(A) is: 

5P(A,)(i?„) = -2 - Re/ / 
'^o V/^o/ Jo Jo 

+Ap{a)Bp[a, (3)*] cos^a saiadadl3. (31) 
According to the Lagrange multiplier rule for inequality constraints (also known as Kuhn- 



Tucker's theorem) 



221 ]. there exists a Lagrange multiplier A > such that, if A is the 



optimum field, we have: 

5F{A) (B) - A 5P{A) (B) = 0, for all B in T^o, (32) 

and 

A [P(A) - Po] = 0. (33) 

In the previous section we have shown that P(A) = Pq, therefore the last equation does not 
give new information. By substituting (130!) and (!3T!) into the Lagrange multiplier rule, we 
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get 



^2 

^0 ^0 



27r 



[Ba{a, P)va + Bpi^oi, /3)f /?] COS a sin a da dp 



— A-^ I — I Re / / [Aa{a, /3) Baia, f3)* + Ap{a)Bi3{a, /3)*] cos^ a sina da dp = 0, 
\ V/^o/ Jo Jo 

for all B in 7^0- (34) 

Because B satisfies f l26|) . it follows that in the first integral we may replace B^ and Bf^ by 
5* and 5^, respectively. Hence, 



Re 



^0 



cos a \/io 



1/2 



cos a \ /io 



1/2 



cos^ a sin adad(3 = 0, for all B in T-Cq. (35) 



This is equivalent to 

"Qmax /"27r 



Jo 



/ \ 1/2 \ / \ 1/2 

_ An ( ) Re(A«) Re(5„) - An I — ) Im(A^) lm(Bo 
cosa \fioJ / V/^o/ 



+ 



COS a 



An(^^^ ^ Re(A^) j Re(5^) - An (^^^ ^ Im(A^) Im(S^) 



X cos^ a sin adadf3 = 0, 



for all B in TYq, i-e. for all B for which 

Qmax /•27r 

/ \hn.{Ba)va + Im(i?^)f^] cos^ a sin adad(3 = 0. 
'0 io 

Choose first 5^ and Bp real. Then (1371) is obviously satisfied and (136!) implies: 

1 / \ 1/2 



(36) 



(37) 



An V ^0 / cos a 
1 / \ 1/2 

An V f / cos a 



(38) 
(39) 

[Im(Ac,) lm{Ba) + Im(Ac,) Im(5„)] cos^ a sin a = 0, (40) 

'0 Jo 

for all Ba, Bp that satisfy (!37|) . This can be stated alternatively by saying that if B = 
Bade + is perpendicular to (wq,/ cos q;)q; + {vp/ cos a)(3, then B is perpendicular to 



By substituting this in (15^ it follows that 
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Im(y4a)a + lm{Ai3)f3 (perpendicular means here of course with respect to scalar product 
We conclude that lm{Aa)6!. + lm{Ap)P is proportional to (f^/ cosa;)Q;+ {v^/ cosa)^: 

lm{A^) = C^, (41) 
cos a 

lm{A^) = (42) 
cos a 



for some constant C. We shall now show that C = 0. By substitution of (|4T|) and (1421) into 
(1261) we get 

C / (f^ + f^) COS a sin a da = 0. (43) 

Jo Jo 

If C 7^ 0, then we must have 

Va = Vi3 = 0, for all a, (3 with < a < ftmax, Q < (3 <2tx. (44) 

Use the expressions (12T|) . fl22l) for i^q, and t;/^ in terms of the cartesian components v^iVy 
and Vz- It is then easily seen that (jSl) implies: Vx = Vy = Vz = 0. This contradicts the 
assumption that v is a unit vector. Hence C = 0. 

We thus conclude that the plane wave amplitudes of the optimum field are given by 

1 / \ 1/2 



^a = T^ - (45) 

An v Co / cos a 

1 / \ 1/2 



Ap = —[^\ (46) 

An \€o J cos a 

The Lagrange multiplier A can be determined by substituting (H5l) and (H6|) into -P(A) = Pq 
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and then using ([22]). We find 

1/2 /-amax /"Ztt 



P(A) 



n / /io 



2A2A2 V eo 



^0 



(f ^ + f ^) sin a da dP 



(47) 



^ r /-"max |-27r 

—-—IT I — I <vl / / (cos^ a cos^ /5 + sin^ /3) sina (ia 
2A2A2 V eo ; I io io 

"Clmax /"27r 

.2 / / /'„^„2 ^, 2 /5 I „^„2 



(cos a sin /5 + cos /5) sin a (iac?/? 



'0 JO 

-v1 / / sin^ a da d(3 
Jo Jo 

sin^ a cos sin /? da 

'0 JO 

/•Qmax r27r 

— 2f ^ / / COS a sin^ a cos (3 da d(3 
Jo Jo 

rOmax /•27r 

—2vyVz / / COS a sin^ a sin (3 da d(3 
Jo Jo 

g cos ftmax - ^ COS^ Omax ) {v^x + 



vm / /io 



2A2A2 Veo 

'4 



g 2 COSttmax + 2 COS ttmax ) 



2A2A2 Veo 



1/2 



(48) 



It follows from P(A) = Po and A > 0, that 



A 



TT ?i 



1/2 



^^o 



1/4 



^ - COS Omax " ^ COS^ Omax " sln^ tt^ax COS Omax^^^ 



1/2 



2pJ/% Veo 

Herewith the derivation of the plane waves amplitudes of the optimum field is complete. 
The maximum of the field component at the origin, i.e. of P, is: 

^2 rOtma.x r^TT 

Fmax = F{A) = / / [Aa{a, /3) Va + Afj{a, f3)vj3] COS a sin a da d(3 
\ Jo Jo 



.(49) 



1 n f fiQ 



AA2 V eo 



'0 

2APo 



1/2 /-amax /•2n 



JO 



(t;^ + v'^) sin a da dP 



■^0 V eo 



4 1 2 .2 2 

— COS ttj-nax g COS Omax SlU Omax COS a^nsLx'^z 



1/2 



(50) 



where we used fl47l) and (l49l). 
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IV. THE OPTIMUM ELECTROMAGNETIC FIELD 



The electric field amplitudes of the plane waves of the optimum field are given by 



= — — [Vc,a + Vfsf3] 

An V Co / V / cos a 



1/2 r 



1 / /ip 

An Veo 
An Veo 



^ ^ sm /3 - \ f . n cos (5 ^ . 

cos pa p f + sin p Q: H p \ Vy — tan a a f ^ 

cos a / V cos a ' 



cos a cos^ (3 + 



sin^ /3 



sin^ a 



+ 

+ 



sin^ a 



cos /? sin /3 X + cos a sin (3 + 



cos a / cos a 
2 ^ , cos^ /5 



cos a 

sin a cos /3x — sin a sin /3y + 



cos a 



cos /5 sin /5 y — sin a cos /5 z 
y — sin a sin /3 z 



sin^ a 



cos a 



(51) 



If we write the right-hand side of (13T]) as the product of a matrix and the vector v on the 
cartesian basis x, y, z, we get 



1/2 



COS a cos^ /3 + 



sin P 



COS /3 sin /3 — sin a cos /5 



cos [3 sin /5 cos a sin^ /3 + 



cos^ /3 



— sin a sin (3 



\ 



sin a cos (3 



— sin a sin (3 



We shall use cylindrical coordinates v?, z for the point of observation r. There holds 



^ = rsm99, 2; = rcosv9, 



(52) 



and the unit vectors (p, g are defined by 



Q = COS X + sin (f y, 
(p = — sin 9? X + cos ip y. 



(53) 



Then 



k ■ r = kon{x sin a cos /? + ?/ sin a sin /3 + z cos a) 

= kon{g cos sin a cos /3 + sin (yj sin asinP + z cos a) 

= A;on[f)sinacos((/? — /5) + zcosa] 

= kongsina cos{(p — P) + konz cosa. 



(54) 
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The optimum electric field in a point r with cylindrical coordinates z is then: 



n 



^0 ^0 



/ \ 1/2 



e* ^ sin a cos a da dfi 



sma e 



ikonz cos a 



f.27r 



1 — sin^ a cos^ /3 — sin^ a cos /5 sin — cos a sin a cos /5 



— sin^ a cos P sin /3 1 — sin^ a sin^ /3 — cos a sin a sin P 
^ — cos a sin a cos /? — cos a sin a sin j3 



sin^ a 



\ 



y^^ikong sin« cos{l3-ip) 



(55) 



Furthermore, the magnetic field amplitudes of the plane waves are (fT9l) : 



1 f Ato 



1/2 r 



sin /5 
cos a 

1/2 



a + cos P(3 ] Vx 



An \^ eo 

cos [3 



1/2 



-VpOL + 



cosa 



cos a 



a. + sin /5/3 — tan a f3vz 



+ [(cos^ P — sin^ /5)x + 2 cos /5 sin /5 y — tan a cos P z] f 



+ (sin /5 X — cos /5y) tan aV;^} (56) 



Hence, on the cartesian basis 



1 A^o 



An \ eo 



1/2 



- A^(a,/3)Q: + A„(a,/3)^ 
— 2cos/3sin/5 cos^ /5 — sin^ /? tana sin /3 
cos^ /3 — sin^ P 2 cos /3 sin /3 — tan a cos P 
— tan a cos P 



tan a sin P 



I 



(57) 



The optimum magnetic field is thus 

-Umax 











-Ap{a^ P)dL + Aa{a, P)(3 e* sin a cos a da 



27r 



v 



— cosa sin(2/5) cosacos(2/3) sin a sin/? 
cosacos(2/5) cosasin(2/?) —sin a cos /? 
sin a sin /? — sin a cos /3 



n 



sinae^^""" ™'"da 



^ikgng sin a cos(/3— yp) 



(58) 
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The integrals over (5 can be computed with the following formulas 23 1 



2-K 



ei(^o<(^~^)cos{mp)dp = 27rrj„(C)cos(mv7), 



2n 



,iCcos(/3-v3) 



sm{mP) dp = 2-71 i"^ JmiC) sm{m(p) , 





for m = 0, 1, 2, . . .. Hence, 

f.27r 



cos^pdp = 71 [Jo(C) - ^2(0 cos(2v?)] 







2n 



□«Ccos(/3-(/3) - 2 
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sin'pdp = n [Jo(C) + -^2(0 cos(2(^)] 
cos /? sin /? = — 7rJ2(C) sin(29?). 



By using the notation 



Jkonz cos a u 



co?y a sin^ a Ji{kong sin a) da, 



the electric and magnetic fields can be expressed on the cartesian basis as: 



n 



E{Q,ip,z) =71- 



1/2 



X < 



I v 










0,3 







2^o' i.Q^z) ) 



' 52'^(£',^) cos(2v9) 5(2'^(^, 2)sin(2v9) -2?5(|'^(^, 2) cosy? ^ 

cf/{Q, z) sin(2v?) -9Vi.Q^ ^) cos(2v9) -2ig];'^{Q, z) sin ip 
^^—2igl''^{g,z)cosip —2igl''^{g,z)sinip J 



Yi{Q,ip,z) 



27in^ 



^ 92^{Q,z)sm{2ip) 
-9Y{.Q^ z) cos(2v9) 

y ^9l [Q,Z)SITILP 



-92^ {q^ z) cos(2v9) 
-9Y{.Q^ z) sin(2v9) 

■ 0,2/ N 

(^,2;)cosv? 



• 0,2/ s . 
• 0,2/ \ 

-ig{ [q,z) cos if 




(59) 
(60) 



(61) 
(62) 
(63) 



(64) 



(65) 



(66) 



Alternative concise expressions are obtained when the electric field, the magnetic field and 
the vector v are all written on the local cylindrical unit basis {q,(^,z} (attached to the 
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point of observation r = gg + zz): 



/ N 1/2 



irn 



\ 1/2 



{ [(s'o'H^',^;) + go^{g,z) + g2^ig,z)) {v^ cos (p + Vy sin cp) - 2igl'^{g, z)v^] g 
+ [9o^iQ,z) + go^ig,z) - g2^ig,z)] {-v^ sin + Vy cos (p) (p 

+ [-2igl''^ {g,z){v^ cos (p + Vy sin Lp) + 2gQ'^{g,z)v,] z} (67) 

{ [(fi'o'^(^', z) + fi-o'^l^i, + g2\Q, z)) - 2igl'^{g, z) t;^] g 
+ [^o'^(^, z) + ^o'^(^, 2) - g^/{g, z)] Cp 

+ \-'^igY'{Qi z) + 2gQ'^{g, z) w^] z} (68) 

H(^i,v2, 2;) = — T-y {fi'2'^(^', 2;) (-t'xsinv^ + t'ycoscy^) ^ 
+ [fl'2'^(^'5 z) {vxcos if + Vysimp) + i5'i'^(^i, z)v^] (p + ig^{^{g,z) {-v^ sinip + VyCosip) z } (69) 
[gl'^io, z)v^g+ [gl'^{g, z) + igi''^{g, z)v^] Cp + ig^{^{g, z)v^z} . (70) 



The time averaged Poynting vector of the optimum field is 
ReS(^,V3,2;) = ^Re[E(^, v?,2;) X H(^,v3, z)*] 



with 



^Re {[E^(^, V?, z)H,{g, (p, z)* - E^{g, (p, z)H^{g, p, z)*] g 

+ [Ez{,Q, P, z)H^{g, p, z)* - E^{g, p, z)Hz{g, p, z)*] (p 

+ [Eg^g, p, z)H^{g, p, z)* - E^{g, p, z)H^{g, p, z)*] z} , (71) 



ReSg{g,p,z) = ^Re[E^{g,p, z)H^{g,p, z)* - E^{g,p, z)H^{g, p, z) 

1/2 

,2 



3 / \ i/^ 

-Im [{gV + ^^o"^ - 9r){9rr] vl + 2Im ^''{9^ 
-2Re [gl^\g'{'r - g'o'igl'')*] v.v^} , (72) 



,0,1 , „2,1 „0,3\/„0,2x*l „.2 I oT™ r„0,3^„0,2x*l ^2 



ReS^{g,p,z) = ^Re[E,{g, p, z)Hg{g,p, z)* - Eg{g,p, z)H^{g,p, z)' 



2 n 



3 / \ 1/2 

, (^^j {Im [-2,}'^(,,^'^)* + (,°'^ + g^ + ^2°'^)(^?'^)1 
-2ReK'^(^?'^)* + ^,°'^(^,^'^)*] v,v^}, 

(73) 
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3 / \ 1/2 

2 n / /io \ r o 0,1 , 2,1 , 0,3\/ 1,1\*1 2 

= ( — 1 |-Re [(^0 +^0 +9i ){92 ) \ 

+Re [{gr + ^7o''^ - 9r){9'/r] < + 2Re 

-Im [(^°'' + 9l'' + ^72'')(^?'')* + ^7!''(^2'')*)] • (74) 

V. OPTIMUM FIELD DISTRIBUTIONS 

We will consider now in more detail subsequently fields obtained by optimizing the lon- 
gitudinal, the transverse and an intermediate component. 

A. Optimum longitudinal component 

In this case: 

V = z. (75) 

Then (l68l) and (1701) become: 

27m ( \ 

E(^, z) = ( ) [-Wi^iQ. z) Q + g^'^ig, z) z] , (76) 

and 

H(^,<^,2) = -^^^^?'^(^,2;)(^. (77) 

Since 5'i'^(0, 0) = 0, the electric field in the origin is parallel to the 2;-axis, hence it is purely 
longitudinal in the origin. In the z = 0-plane the functions g'^''^ are real and therefore (176|) 
implies that the polarization ellipse of the electric field in that plane has minor and major 
axis parallel to the Q- and z-axis (which is the major and which the minor axis depends 
on the relative values of go'^{g,0) and 5'i'^(f?, 0)). In the z = 0-plane, the phase of E^ 
is ±7r whereas the other electric field components have phase ±7r/2. The magnetic field 
is everywhere parallel to c^, i.e. it is azimuthal. In Fig. [1] the normalized distribution 
\Ez{x,y ,0)1"^ of the optimum field with maximum longitudinal component for NA/n = 0.5 
in the z = 0-plane, is compared to the normalized total electric energy density: 

\E,{x, y, 0) |2 + \Ey{x, y, 0) ^ + \E,{x, y, (78) 



17 




-1.5 -1 -0.5 0.5 1 1.5 ■ -1.5 -1 -0.5 0.5 1 1.5 



Figure 1: Left: the normalized distribution of \Ez\'^ in the z = 0-plane for the field with maximum 
longitudinal component when NA/n = 0.5 and Pq = I {X = Xo/n is the wavelength in the material 
with refractive index n). For comparison, the normalized electric energy distribution (I78p of the 
j;-polarized focused plane wave is shown at the right for the same numerical aperture. The total 
flow of power in the z-direction is the same for both fields. 



in the focal plane of the focused x-polarized plane wave. The coordinates x, y are expressed 
in units of the wavelength A = Ao/n in the material with refractive index n. The Formulas 
for the electric field of a focused plane wave are given in the appendix (see IA14p . flAlSp . 
(1A16P ). The power flow in the ^-direction of the optimum longitudinal field and the focused 
plane wave are the same. The distribution of the optimum \Ez\'^ is rotationally symmetric 
while that of the electric energy density Eq. fl78l) of the focused linearly polarized plane wave 
is elliptical, with short axis parallel to the ^/-direction (i.e. perpendicular to the direction 
of polarization of the focused plane wave). In Fig. [2] the distributions are compared for 
NA/n = 0.9. In this case the short axis of the elliptic distribution is considerably shorter 
than the long axis. In Fig. [3] the squared modulus \Eg\'^ and |Ep are shown. Due to the 
broad doughnut shaped distribution of the radial component, |Ep is broader than the Airy 
spot. 

Cross-sections along the short and long axes are shown in Fig. IHfor both NA/ra = 0.5 
and NA/n = 0.9. To compare the shapes the maxima of all cross-sections are rescaled to 
1. It is seen that the longitudinal component has smaller FWHM (Full Width at Half 
Maximum) but also higher secondary maxima. 

The FWHM of the optimum longitudinal component is for NA/n = 1 almost identical to 
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Figure 2: Same as Fig. [T] except that now NA/n = 0.9. 
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Figure 3: The squared modulus \Eg\'^ of the radial component and electric energy density |Ep of 
the electric field with maximum longitudinal component for NA/?i = 0.9. The maximum values 
of \Eg\^ is approximately 25% of the maximum of the squared modululus of the longitudinal 
component. 



that of the longitudinal component in [18[, obtained by focusing a radially polarized beam 
using a ring mask function (with radius 90% of the total pupil). However, the side lobes are 
higher at the cost of the central maximum compared to our longitudinal component. This 
is of course not surprising because the longitudinal component in was not optimized for 
a high maximum on the optical axis. 

In Fig. [5] the FWHM of ?/, 0)p of the optimum field is compared to the FWHM 

in the x- and y-directions of the electric energy density |E(x, y, z = 0)p of the focal spot of 
the x-polarized plane wave. At the left, the FWHM is shown as a function of the numerical 
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Figure 4: Cross-section of the rotational symmetric \Ez\^ (solid curves) of the field with maximum 
longitudinal component and cross- sections of |E(x,0,0)|2 (dashed) and |E(0,y,0)p (dotted) of the 
Airy spot for the x-polarized focused plane wave, for NA/n = 0.5 (left) and NA/n = 0.9 (right). 
\ = \o/n and the maxima of all quantities are rescaled to unity. 

aperture NA/n, at the right as a function of n/NA, both in units of wavelength. It is 
seen that the FWHM of the longitudinal component is smaller than the FWHM in both 
the X- and y-directions of the focused spot. The FWHM as a function of n/NA is almost 
linear, but the slope is smaller for the longitudinal component. This means that the spot 
size of the optimum longitudinal is relatively smaller than the Airy spot when the 

numerical aperture is larger. Nevertheless, for all values of the numerical aperture, the 
longitudinal spot is narrower than the Airy spot. As shown in Fig. [6l the maximum 
amplitude |i?2(0)| of the optimum longitudinal component is for most values of the numerical 
aperture smaller than the maximum amplitude \Ex{0) \ of the focused x-polarized plane wave. 
But for NA/n > 0.65 this ratio is already more than 0.5 and for NA/n > 0.994 the maximum 
longitudinal component is even larger than the maximum amplitude |-E2;(0)| of the focused 



20 




Figure 5: FWHM of \Ez{x,y,0)\'^ of the field with maximum longitudinal component and FWHM 
of |E(x,?/, 0)p along the x- and y-directions of the focused x-polarized plane wave. At the left the 
FWHM is shown as function of NA/n, at the right as function of n/NA in units of A = Ao/ra. 




NA/n 

Figure 6: Ratio of maximum longitudinal amplitude \Ez\ and the amplitude \Ex\ of the focused 
x-polarized plane wave for the same power and as function of NA/n. 

x-polarized plane wave. 
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Along the optical axis we have ^ = and 

/"Qmax 

^0'3(o,^) = / e^'^'^^^^'^sin^arfa 







i ( 2 



cos a gifco"'^ cos Ctmax ^ikQUZ 

COS^ ^ikonz COS a^iz.^ ^ikonz 

k^nz 



(79) 



Because Ji(0) = 0: 



g]^\Q^z)=g\:\Q,z) = Q. (80) 

It thus follows from (176!) that along the optical axis the electric field is parallel to the ^-axis 
and that the modulus of the field is symmetric with respect to the focal plane. Furthermore, 



|E,(0,0,0)P \gl^\Q 



(81) 



where 



2 1 

5(0'^ (0, 0) = - - cos ttmax + - COS^ a^ax- (§2) 



We define the focal depth of the optimum longitudinal component as the distance Az to the 
focal plane for which the ratio 

^ 0.8. (83) 

For the focused linearly polarized plane wave, the electric energy density on the optical 
axis is given by flA40p : 

|E(x = 0,y = 0,2;)p = \E,{x = 0,1/ = 0,2;)^ + \Ey{x = 0,|/ = 0,2;)^ + \E,{x = 0,y = O,^;)^ 

2 2 -f 2 

= ^^l4''(0,.)+^|'\0,^)r. (84) 

^0 



In the origin we have (1A41I) 



|E(0,0,0)p 



^0 



2 2 1 ^ 

- (1 - cos^/^ a„iax) + 5(1- cos^/^ ama.) . (85) 



The energy density on the optical axis is again symmetric around the focal plane and the 
focal depth is the distance Aza such that 

|E(a; = 0,y = 0,A^^)|2 



|E(0,0,0)|2 
22 



0.8. (86) 
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Figure 7: Focal depth in unit of A = Ao/^^ as function of NA/n for the optimum longitudinal com- 
ponent and for the energy density of the focused linear polarized plane wave. The approximative 
paraxial focal depth: O-bXan/NA^ is also shown. 

Both Az and Aza are shown as functions of NA/n in Fig. [71 The approximative focal 
depth 0.5Ao?T^/NA^, which is vaUd in the scalar paraxial theory, is also shown. It is seen that 
for NA/n < 0.4 the focal depths of the focused linearly polarized plane wave calculated in 
the vectorial and the scalar paraxial theory are almost identical. The focal depth of the 
optimized longitudinal spot is for NA/n < 0.8 larger than that of the Airy spot, while for 
NA/n > 0.8 it is smaller. 

Eqs. ([72]), dZl) and dUD imply 

2^2^3 / \ 1/2 

S(^, (fi, z) = l^-^j {im [go^ig, z){g'l'^{g, z))*] g + Re [gl''^ig, z){g'^'^{g, z))*] z }(87) 

The yj-component of the Poynting vector thus vanishes and the Poynting vector is indepen- 
dent of the angle ip. In the {z = 0-)plane the energy flows in the 2;-direction. In Fig. [8] the 
time averaged Poynting vector of the field with maximum longitudinal component is shown 
in the {z,g)-p\a.ne for the cases NA/n = 0.5 and NA/n = 0.9. The (/^-component of the 
Poynting vector vanishes and the other components and are independent of rotation 
angle ip. The Poynting vector vanishes on the 2;-axis. 
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zlX zlX 

Figure 8: Poynting vector S = SqQ+Sz^ for the fields with maximum longitudinal component, 
in the (z, x)-plane (i.e. the (p = 7r/2-plane) for NA/n = 0.5 (left) and for NA/n = 0.9 (right). The 
Poynting vector is independent of the angular coordinate (p. 



B. Optimum transverse component 

We choose the transverse component parallel to the x-axis: 

V = X, (88) 

i.e. the x-component of the electric field in the origin is optimized. On the local basis 
(attached to the point of observation) of cylindrical coordinates v?, z we have: Vg = cos if, 
Vy, = — simp, Vz = 0. Hence (165!) . (!68|) . (!66|) and (ffOj) imply that the optimum electromag- 
netic field becomes: 

E(^, (f, z) = i-^j { [go^ig, z) + g^^^Q, z) + g^/[Q, z) cos(2(^)] x 

+g2^i0, z) sin(2v9) y - 2ig^{'^{Q, z) cos (/? z} (89) 
/ \ 1/2 

nn fiQ \ .r 0,1/ X 2,1/ \ I 0,3/ \-\ 

{[9o {Q^z)+gQ' {g,z)+g2 {g,z)\ cos ipg 



- [go^iQ,z) + gl'^{g,z) - g^'^ig^z)] sincp (p - 2igl''^{g, z) coscpz} , (90) 



ii{g,'^,z) = --TTT [al'^g, z) sm2<p> ± - gl'\g,z)cos2ip y -igi'^{g,z) simp z] (91) 
= TT2" [dl'^i^^^) ^^"^^Q- al'^iQ^z) cos if (fi + ig^''^{g,z) siny^z] . (92) 
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Figure 10: Normalized distribution of \Ez\'^ and of the electric energy density |Ep for the field with 
maximum component NA/n = 0.9. 

The squared amplitudes of the electric field components and the electric energy density in 
the z = 0-plane are shown in Figs. [9] and [TO] when NA/n = 0.9. These field components are 
very similar to the components of the focused x-polarized plane wave for the same NA/n. 
In Fig. [TT] cross-sections along the x- and y-axes of l-ExP for the maximum transverse field 
are the focused linearly polarized plane wave pattern are compared for NA/n = 0.9. The 
maxima of all functions are normalized to 1. Along the y-axis, the E^-spot of the optimized 
transverse field is slightly narrower than for the Airy spot, but has also higher secondary 
maxima. It can be shown that, just as for the focused x-polarized plane wave, of the 
optimized field is widest along the x-axis. In Fig. [12] the FWHM as function of NA/n of 
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Figure 11: Cross-sections of \Ex^ along the x— (left) and ?/— axis (right) for the field with maximum 
x-component and for the focused x-polarized plane wave with the same power. The numerical 
aperture is NA/n=0.9. The maxima of all cross-sections are rescaled to 1. 

\Ej}^ of the optimized field is compared to the FWHM of |i?xp of the x-polarized plane 
wave. The field distributions have elliptical shape. The FWHM is therefore defined with 
respect to the radial distribution obtained by averaging \Er^{^Q^if^z = 0)p over the angle 
< < 27r, as explained in flA26P in the Appendix. Note that, in contrast to Fig. [5l 
where the optimum component was compared with the intensity of the Airy spot, we 
here compare the squared moduli of the x-components of the electric fields. It is seen that 
the FWHM of the optimized and the focused plane wave are almost identical. 

The optimized electric field at the focal point is pointing in the x-direction and is given 

by 

(93) 



E,(0,0,0) = 



cos Omax g COS Omax 



■^0 V ^0 

Formulas for the field distributions of the focused x-polarized plane wave are derived in the 
Appendix. For a power given by flA30p . the x-component of the electric field in focus is 
given by (1A41I) . Hence if the power is Pq "we have for the focused x-polarized plane wave: 



^0 V Co / 



1/4 



1 /16 2 

cos^' ar 

smamax V15 3 



-cos5/2 



Or 



(94) 



The ratio of 



and 



( 3 COS ftmax 3 COS Omax) 
Sm CI max ^[g 2 
15 ~ 3 



1/2 



^ - ^ COs3/2 ameo, - I COs5/2 



(95) 
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Figure 12: The FWHM of along the X-axis of the optimized transverse field and the FWHM 
of the averaged of the x-polarized focused plane wave. 

is shown in Fig. [T3]as function of NA/n = sinamax- For NA — ^ the ratio becomes 1, while 
for NA/n — > 1 it becomes 15/(8v^) = 1.0825. Hence, the optimized field can have more 
than 8% higher amplitude than the x-polarized Airy spot of the same power. 

C. Intermediate case 

Without restricting the generality, we may assume that the vector v is in the (x, z)-plane, 
i.e. Vy = 0. It is then convenient to write the optimum focused electric and magnetic field 
on the orthonormal basis with unit vectors v, y and v x y. We have 

V = + v^z, V X y = -w^x + t;^z, (96) 

and therefore 

X = t;x.v - fz(v X y), z = v^ir + v^(y x y). (97) 
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Figure 13: Ratio of the amplitudes of Ex{0,0,0) at the focal point of the field with maximum 
|£^x(0, 0,0)1 and of the focused x-polarized plane wave with the same power. The ratio is for 
NA/n = 1 maximum with value 15/(8\/3). 

Hence, (l65l) and (l66ll become: 

^ / \ 1/2 

E(^, ^, z) = TTjj^ l^-^j { [(^o'^ + ^o'^ + ^2'^ cos(2(/9)) Vx - 2igl'^ cos ip v,] x 

+ [g2^ sin(2v?) - 2ig\''^ sin ipv^]y + [-2ig\'^ cos (p + 2gl'^ v^] z} 
/ \ 1/2 

/ ^0 \ r r/ 0,1 , 2,1 , 0,3 /o \N 2 /I • 1,2 ID 0,3 21 - 

= (^~J iK^o +^0 +92 cos{2ip})Vx- Atg{ coscpv^v^ + 2gQ' v^\v 

+ [5(3'^ sin(2(y9) i;^ - 2igl''^ sin f^] y 

I r o- 1.2 2 / 0,1 , 2,1 o 0,3 , 0,3 /o ^^ I o- 1.2 21 - 

+ [-2z5(i' cos (y? t;^. - (5(0' +5-0 -25-0 +92 cos{2ip)) v-^v^ + 2ig{ cosipv^\ v 



and 



H{g,ip, z) 



2'Kn^ 
AAl" 



{ [5(2'^ sin(2v9) + igi"^ sin (yj x 



-[5(2' cos(2v9) t;^. + cosc/^i;^] y + zgfi' sin^jw^^z} 
^TTT { [^2'^ sin(2y?) + 2ig^{^ sin f J v 

AAq 

- [^fg'"^ cos(2v9) t;^ + igl'^ cos </? f ^] y 

I r ■ 0,2 . 2 1,1 • /o N ■ 0.2 . 2i - , - 1 



(99) 
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Because gi''^{0, z) = when £ > 1, it follows that on the optical axis: ^ = 0, we have 

E(0, yp, z) = vr^ { [{g'/iO, z) + ^^''^(O, z)) vl + 2^°''(0, ^) ^.^] v 

- [^7|;''(0, ^) + ^7o''(0, ^) - 2^?|;''(0, z)] v^v, V X y} , (100) 

while the magnetic field vanishes there. It follows in particular that the electric field in the 
origin = z = is parallel to v only when v is perpendicular to the optical axis (f^ = 0, 
transverse case) or parallel to it iv^ = longitudinal case). In all other cases the projection 
of the electric field at the origin on the plane perpendicular to v is nonzero. 
The u-component of the electric field is 

^ / \ 1/2 

B{g,(p,z)-\r = TT——(-^\ {[g^'^ + gl'^ +g^'^cos{2ip)] vl-Aigl'^ COS if v^v, + 2gQ'^vl} . 

(101) 

In particular, since the functions g'^''^ are all real for z = 0, we find for the squared modulus 
of the w-component in the z = 0-plane: 

iTji/ n^ ' |2 2 f f^o\ fu 0,1 I 2,1 0,3\ 2 , o 0.3 2i2 

|E(^,(^,0) ■ v| =^J^{—){ly9o +9o - 92 )v^ + 2go' v^l 

I A [/■ 0,1 , 2,1 0,3\ 0.3 4 I o/ 0>3 0,3 , r,/ 1>2\2\ 2 21 2 

+4 [(^0 +^0 - 92 )92 ^x + 2(^o 92 +2(^i' ) Kv^\ cos if 
+%2°'')'^^osV}, (102) 

where all gi''^ are evaluated at g and z = 0. If 'dv is the angle between the vector v and the 
positive 2-axis, we have 

Vx = sin'dy, Vz = cos'dv. (103) 

In Figs. UM and [T5| |E ■ vp, {Eyl"^, |E ■ (v x y)p and the electric energy density |Ep are 
shown in the z = 0-plane for the optimum field with = 30° and NA/n = 0.9. Figs. [16] 
and [17] correspond to 'dy = 60°. 

In Table 1 the maxima of the indicated field components in the focal plane are listed for 
the optimum fields corresponding to NA/n = 0.9 for several choices of 'dy. In all cases the 
total fiow of power in the pupil of the lens is 1 W. 

Eqn. (11021) is a quadratic function of cos^ (f and since the coefficient of cos^ (f is non- 
negative, it follows that, for every g, the maximum is attained when cos^ ip = 1, i.e. when 
if = and when if = n. To obtain a useful measure of the spot size, we average (11020 over 
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< < 27r: 

'•27r 

27r 



1 f iTrf -|2 J 2 f /^oA f-r/ 0,1 , 2,1 0.3\ 2 , o 0,3 2i2 

- / |E(^, if, 0) ■ v| = IT [^) I t(^o + ^0 ~ ^2 + 2^0 

+ l(9rfv'\. (104) 



where we used 

1 /"^'^ 1 1 /"^'^ 3 

— / cos^ (iy? = - , — / cos"^ = -. (105) 

2ti Jo 2 27r Jo 8 

The FWHM in the z = 0-plane is then defined as the value 2gQ such that 

r2TT 

In Figs, and the FWHM of and the maximum |£'^(0)| are shown as as function of 
for several values of the numerical aperture and for Po = 1 W. In Fig. the amplitude of 
I £"2; (0)1 of the focused x-polarized plane wave is shown for the same power are indicated by 
dots. In contrast to the FWHM, the maximum Ey{Q) is a monotonically increasing function 
of 



^ [ \E{go, 0) ■ v|2 dif = l |E(0, 0, 0) ■ v| 
271 Jq 2 



(106) 



VI. THE OPTIMUM FIELD IN THE LENS PUPIL 

The optimum fields in the z = 0-plane can be obtained by focusing appropriate pupil 
distributions. We derive these pupil distributions in this section, using the vector diffraction 
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x/X x/X 

Figure 16: The normalized distribution of \Ey\'^ (left) and \Ey\'^ (right) in the z = 0-plane for the 
field with optimum field with i?„ = 60° when NA/n = 0.9. 



theory of Ignatowksy and Richards and Wolf. This theory is based on 1. Debye's approxi- 
mation which expresses the plane wave amplitudes in image space to the field vectors in the 
entrance pupil, and 2. the Abbe's sine condition to guarantee conservation of energy. 

Consider a beam that is incident on the pupil of a diffraction limited lens with numerical 
aperture NA. The electric field in the pupil is given by 

^"{xp, Vp) = EP{xp, yp)<pp + E^ixp, yp)Qp, (107) 

where Xp and t/p are cartesian coordinates in the lens pupil, which are parallel to the x- and 
y-coordinates of the cartesian system {x, y, z) in the focal region, and ipp and Qp are polar 
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Figure 17: Normalized distribution of |i? • (v x y)P and of the electric energy density |Ep of the 
field with maximum u-component for = 60° and NA/n = 0.9. 



[degrees] 


l^.l [V/m] 


\Ey\ [V/m] 


|^-(vxy)| [V/m] 


90° 


45.28 


6.03 


13.67 


60° 


42.93 


7.33 


13.81 


30° 


37.76 


14.17 


14.04 


0° 


34.87 


17.72 


17.72 



Table I: Maxmima of the optimum electric field components (not of their squares) in the focal 
plane when NA/n = 0.9 and for several choices of the angle The total power in the lens pupil 
is in all cases unity. Note that the maximum of \Ey\ and |i? • (v x y)| are not attained at the focal 
point. 



coordinates in the pupil: 



Vp = Qp sin (pp, 



(108) 
(109) 



with unit vectors 



Qp = cos(y9pX + sinv?py, 
iPp = - sin LppX + cos ippj. 



According to the theory of Ignatowsky 



and Richards and Wolf 



a, a 



(110) 

(111) 

the field in 
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Figure 18: FWHM expressed in units of A = Ao/?t. of the optimum v-component as function of i}.. 
for several values of NA/n. 

50 
45 
40 

, — . 35 
> 30 
2, 25 
^^'^20 

15 
10 

^0 10 20 30 40 50 60 70 80 90 100 

e 

V 

Figure 19: E^{0) in V/m of the field with optimum v— component as function of the angle between 
V and the z-axis (i.e. v = sini?„x + cosi?^z), for several NA. The power incident on the lens is 
Pq = I W. The dots are the amplitudes |£'a;(0)| of the focused x-polarized plane wave for the same 
power and the same values of the numerical aperture. 
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the pupil of the lens and the amplitude of the plane waves in the focal plane are related by 



where / is the focal distance of the lens and ^(k) is a matrix which rotates the electric field 
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in the pupil plane in the direction perpendicular to the wave vector k = kon{sm a cos (3 "k + 
sin a sin P a + cos a /3) : 



^{k)EPixp,yp) = -EP{xp,yp)ci-EP{xp,yp)(3, (113) 

with 

k 

Xp = Qpcosipp = -f-^, (114) 

k 

Vp = Qp sin ipp = -f-j^, (115) 

and a and f3 are given by ([8]), (Q. Since a and /3 are the polar and azimuthal angles of the 
wave vector k: 

Qp = /sin a, (116) 

iPp = P + 'K, (117) 

Hence, expressed in terms of a and /?: 

/ 

A(a,/9) = —2ni- TZja, l3)Ei^ (— f sin a cos 0, — f sin a sin /j) 

Konvcosa — 

= —2ni- r_E^ (— / sin a cos/5, — / sin a sin/3) a. 

koTiy/cosa ^ 

+E^{—f sin a cos P,—f sin a sin P) (3 . (US) 

Conversely, given the plane waves amplitudes A{kx, ky) = AaOt + ^/3/3 of the field in the 
focal plane 2; = 0, the corresponding field in the pupil of the lens is: 

E^(^P, Vp) = A (^o^)'^'^'^V (k)-iA(fc., ky), (119) 
27r / — 



with 



kx = -fco?^y, ky = -koUj-, (120) 



(x2 + y2)i/2 < /NA/n and 

^(k)-'A = -A,^p-A^c^p. (121) 

In terms of a, (3 we get: 

w{qp,^p) = ±^ cos'/' an{<y,Pr'M»,P), (122) 
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where Qp, and a, (3 are related by flllGp and (11171) . By substituting the expressions for 
the optimum plane wave coordinates (HSj) . fH6l) : 



we find 



where 
A 



Ad 



1 //LtO 



1/2 



An \ eo / cos a 

1 / \ 1/2 
1 / /iO \ Vfi 



An V eo / cos a 



2 /Cgn 



1 //.o\^/^ 1 



AAo/ V^o/ a/cos a 



TT n 



1/2 



2 pV2), 



(123) 
(124) 



(125) 



cos Omax - g COS^ ttmax " slu^ O^ax COS Omaxf ^ 



1/2 



(126) 



Using (EH), (El and 



cos a 



/ 



sm a 



Op 



(3 = <P: 



(127) 

(128) 
(129) 



we obtain: 



AAo/ V eo / 



\ 1/2 

-) [("■ 



COS a cos P + Vy cos a sin — t;^ sin a) q„ 



+ {-v^ sin /3 + t;^ cos (3) (fip] 



1 



a/ cos a 



\ 1/2 

AAo/ V eo / 



f 



cos (pp + Vy sin ipp) + 



Jl/2(y2 _ ^2)1/4 



Op 



sin + Vy cos V9p) (^^ 



AAo/ V eo 



1/2 



/ 



^p 



J1/2(J2_ ^2)1/4 



^p + 



^1/2 



"Pp 



(130) 
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where u^, are the components of the vector v on the polar basis Qp, (p^ in the pupiL It 
is seen that the optimum electric field in the pupil is linearly polarized with direction of 
polarization and amplitude that depend on the radial coordinate only. 
In the case of the maximum longitudinal component: Vg = 

1/2 



0, 



AAo/veo; mp 

whereas for the maximum x-component: = cos(/?p, = 

1/2 



E^(^p, ^p) 



AAo/ V eo 



f 



COS (fp Qp 



AAo/ V eo 



1/2 



— cos (fp + 



- sin Lpp,v^ = 
iP - QlY" 

fl/2 



1, this becomes 
(131) 




0: 



sm ifp 



sm ifp X 



(132) 



One may have expected that the maximum longitudinal component is obtained by con- 
centrating the radially polarized pupil field in a small annular ring at the rim of the pupil. 
Instead, the optimum field amplitude (11311) increases continuously from zero at the center to 
a maximum at the rim of the pupil. This is due to the fact that the total power is an integral 
over the pupil of the squared amplitude of the field while the longitudinal component in the 
focal point is the integral over the pupil of the field itself. If the radially polarized pupil 
field is concentrated in a ring of width 6gp = /NA 6a, the longitudinal field component in 
the focal plane becomes approximately 



Ez{g,(f),0) oc Jo{kongNA) {6a 



,1/2 



(133) 



The longitudinal component thus becomes in the limit 6a ~* proportional to jQ{kongNA) 
and hence decreases slowly as l/V^of? for increasing distance g to the optical axis. To 
keep the total power finite, the amplitude vanishes with the square root of the width of the 
annular ring. The FWHM of the squared amplitude decreases with the width of the ring 
and has minimum value of 0.35Ao/(riNA) for 6a 0. But for small widths the side lobes 
are relatively strong, the first maximum adjacent to the center having squared amplitude 
that is 16% of the central maximum. 

Suppose now that the vector v is in the (x, 2;)-plane and the angle with the positive ^-axis 
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is 'dv Then f llOSp holds and fll30p becomes: 




cos ifp sin + 



/1/2(J2 _ ^2)1/4 



cos-i?. 



V 



(134) 



p 



The pupil field is always linearly polarized, with direction of polarization that varies with 
the pupil point. The phase of the field is the same in all points of the pupil. In Fig. [20] we 
show snapshots of the electric field in the pupil for §^ = {] (i.e. the longitudinal component 
in focus is maximum), = 90° (i.e. the x-component in focus is maximum) and the 
intermediate cases ■i?^ = 30° and = 60°. In all cases NA/n = 0.9. For i)^ = the electric 
field is radially polarized, for -t^y = 90° it is the field of a plane wave polarized parallel to 
the X-axis. For -dy = 30° and i)y = 60° intermediate cases occur. In Fig. [21] the squared 
modulus of the electric pupil field is shown as function of the normalized pupil coordinate 
Qp = Qp/O'i with a = /NA/n, for several values of NA/n and for -^y = 0° (i.e. the pupil fields 
correspond to fields with maximum longitudinal component). The square modulus shown 
has been rescaled by the factor a^, hence the function 



is shown. The reason is that the power in the pupil is proportional to the following integral 
over the normalized radial pupil coordinate: 



It is seen that for higher NA/n the electric energy density is more concentrated near the 
rim of the pupil. But even for NA/n = 0.9, the energy density is smoothly varying inside 
the pupil. Hence it seems quite possible to realize these pupil fields using appropriately 
programmed spatial light modulators. It may be important in some applications to realize 
the radially increasing amplitude without absorbing a substantial part of the light. This 
could be done with optical elements which refract energy close to the axis towards the rim 
of the pupil. 



(135) 




(136) 



VII. DISCUSSION 



We have studied fields which, with respect to some cartesian coordinate system {x,y,z), 
consist of plane waves that all propagate in the positive 2;-direction in a homogeneous 
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Figure 20: Snapshot of the electric field in the pupil that, when focused, yields the field with 
maximum u-component in focus when NA/n = 0.9. The upper left figure corresponds to "Q^ = 0°, 
the upper right to 19^, = 30°, the lower left to ■i?^ = 60° and the lower right to "d^, = 90°. 



isotropic lossless medium with refractive index n. The sinus of the maximum angle be- 
tween the plane wave vectors and the 2;-axis is the numerical aperture of the plane wave 
superposition. Let v be a unit vector. For given v, given NA and given mean power flow in 
the positive ^-direction, we have determined closed formulas for the plane wave amplitudes 
for which the projection of the electric fleld in the direction of v in the origin is larger or 
equal than that of all flelds with the same NA and the same maximum power flow. Closed 
formulae have been derived for the plane wave amplitudes. By choosing v parallel to the 
z-axis, we get the fleld for which the amplitude of the longitudinal component is maximum. 
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Figure 21: Electric energy density a^|E(^)j,p in the lens pupil as function of the radial coordinate 
Qp, corresponding to the maximum longitudinal component in focus for several NA/n. All values 
are relative to the value for NA/n = 0.9 at the rim of the pupil. 

The width of this longitudinal component is considerably smaller than that of the Airy spot 
for the same power and the same NA. By choosing v perpendicular to the 2;-axis, along the 
X-axis say, we obtain the field for which the amplitude of the x-component of the electric 
field in the origin is maximum. This field is similar to the field of a focused linearly polarized 
plane, but it is not identical to it. Also for directions of v that are intermediate between 
longitudinal and transverse, the optimum fields were studied. In particular the FWHM of 
the components in the plane z = where compared for several values of the numerical 
aperture. 

By using the vectorial diffraction theory of Ignatowsky and Richards and Wolf, the field 
distributions in the pupil of a diffraction limited lens of given NA was derived such that the 
focused are the previously discussed fields for which the projection of the electric field along 
the v-direction is maximum in the focal point. Closed formulas were obtained. When v is 
chosen equal to the unit vector along the z axis, the pupil field is radially polarized. When 
V is parallel to the x-direction, the pupil field is almost identical to that of a plane wave that 
is polarized parallel to the x-direction. In general, the pupil fields that yield the optimum 
fields near the focal plane, are linearly polarized and in phase in all points of the pupil. But 
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the state of polarization and the amphtude varies with position in the pupiL 



Appendix A: THE FOCUSED FIELD OF A LINEARLY POLARIZED PLANE 
WAVE 

We recall here the expression for the electric field near the focal plane of a lineary polarized 
plane wave in the pupil of a lens of numerical aperture NA, as derived in the theory of 
Ignatowsky and Richards and Wolf [3|, [J]. With respect the polar basis, the x- 

polarized pupil field can be written 

E^(^p, ^p) = E^iQp. ^p)Qp + EUqp, ipp)(fp, (Al) 



with 



Qp = COS X + sin y, (A2) 
<Pp = - sin X + cos y, (A3) 



and 



E^g{Qp,^p) = COS ifp, (A4) 
E^{gp,(Pp) = -sinv^p. (A5) 

According to (fT8l) . (fT9|) . the plane wave expansion of the electromagnetic field near the focal 
plane can be written as 

,2 /-Ctmax /•27r 



^2 /•Ctmax r^TT 

E(r) = -72/ / {AaCt + A/s^) sin a cos a e'^ '' da d(3, (A6) 
■^0 Jo Jo 

/ e \ /""max r2Tr 

H(r) = -2 — / / {-Af3dL + Aa^) sin a cos a e^^ '' da dp, (AT) 
Aq V/^o/ Jo Jo 

where a and (3 are spherical coordinates for the wave vectors of the plane waves: 

k = kou (cos P sin ax + sin /5 sin ay + cos az) , (A8) 

and d and /3 are given by ([8]), IQ. With respect to cylindrical coordinates: r = gQ + zz, we 
have 

k ■ r = kon [g sin a cos(/3 — ip) + z cos a] . (A9) 
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By substituting dM]), (jM]) into ffTTSD and using ffTTTD : = /9 + vr, we find 

/ cos P 



Ap{a,(3) = -2Txi 



kou -y/cos a ' 
/ sin/3 



/co?T- a/cos a 



(AlO) 
(All) 



Hence, 

A(a,/3) = 27ri 
= 27ri 

and 



f 



kQn^y cos a 
f 



cos P a. — sin P f3 

{ (cos a cos^ /5 + sin^ /?) x + (cos a — 1) cos P sin /?y — sin a cos z} , 

(A12) 



27ri 
27r2 



/ 



koUy/cosa . 
f 

koTiy/cosa 



sin Pa. + cos /9/3 

[(cos a — 1) cos P sin /5x + (cos a sin^ /5 + cos^ P)y — sin a cos /3z] 



(A13) 



where / is the focal length. Substitution into plane wave expansion f lA6l) and using formulas 
§MM and dSH) yields 
zvrn/ 



Eyir) 



EJr) 



■ 1 1 3 ^ 

fl-l' (Q^z) -92 {Q^^) 



and 



i^x(r) 



Hy{r) 



HJr) 



Ao 
innf 

2?™/ 1,2, - 

— — {g,z) cos if, 

^0 



gS' (9, z) + {g, z) + {g, z) - g^' {g, z) cos2v9L (A14) 



sin2<y9, 



(A15) 
(A16) 



ivrnV |eo" f i i 



— [^1' {Q^z)-g^' (^, 2;)]sin(2v?), 



Ao V A^o 

27™V /i;^ 1 2 



il i I 1 I ±1 ^ 

9o' {q, z) + go {g, z) - [g^' {g, z) - g^' {g, z)] cos(2(/?)| 



Ao V A'o 



— g{ [g,z)smip. 



(A17) 

(A18) 
(A19) 
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The cylindrical field components are 

Eg{r) = Ex{r) cos (p + Ey{r) sirup 



and 



■ 1 1 3 1 1 1 3 1 

9o' {q, z) + g^' {q, z) + {q, z) - g^' {g, z) 



innf 

~Ex{r) sirup + Ey{r) cosip 



coscp, 



iirnf 



I I 3 ^ 1 ^ 3 

9o' (Q^ + 9o {q, z) - ^1' {q, z) + g^' {g, z) 



Sm <y9, 



(A20) 
(A21) 



Hg{v) = H^{t) COS (p + Hy{T)sm.ip 

%7TTl^ f I 6q i-l ^1 — \ A \ 

= \ — 9o' {Q,z)+g§' {g,z)+g^' {g,z)-g^' {g, z) simp, (A22) 

Aq \ flQ I J 

H^{r) = — Hx{r) sinip + Hy (r) cos ip 
_ inn'^f 17^ 

In the plane 2; = all functions g'^'^ are real. Therefore, the squared modulus of the electric 
field in the z = 0-plane is 



9o^ (q, z) + g^' (g, z) - g^' (g, z) + g^' {g, z) cosip. (A23) 



|E(^, ip, 



^0 



9f{Q, 0) + gf{g, 0) + gf{g, 0) - gfig, 0) 



1 2 



cos^ ip 



+ 
+ 



9riQ, 0) + grig, 0) - grig, 0) + gfig, 0) 
2 



sin^ (yj 



cos^ ip 



(A24) 



The spot shape is elliptical. To obtain a measure of the spot size we average the squared 
modulus over < ip < 2n: 



2n 



— I \E{g,ip,0)\'dip 



2X1 



9fiQ, 0) + 9 fig, 0) + gfig, 0) - gfig, 0) 



+ 



+ 



9T iQ,0)+9~r iQ,0) -9r iQ,0)+9r (f?,0) 



9r iQ,0) 



^0 



9l'\Q,0)+gl'\g,0) 



+ 



9f{Q,0)-gfig, 0) 



2 1 
+ 2 



9!' (f?,0) 



The FWHM in the z = 0-plane is then defined by 2 go with go such that 

- \Eigo,^,0)\'d^ = -\Em'- 
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(A25) 



(A26) 



The real part of the complex Poynting vector in cylindrical coordinates is: 



Re5,(r) = hie[E^H: - E^H;] 



^0 




Im 



(A27) 



Re5^(r) = ^Re [E^H; - E,H:] = 0, 



(A28) 



1 



Re5,(r) = -Re[E,H;-E^H;] 



2 

2Aq Y /.io 



1.1,2 



gi 



(A29) 



The total power flow P through a plane z = constant is most easily calculated by substituting 
flMOj) and fIXTT]) into One flnds 



2 




sm a„ 



(A30) 



We have along the 2;-axis: 



for all u and and 



,tecosa cos^/^asinac/a 



^ikonz gifcon2 cos Omax 



ik^nz 



COS ^ ttmax 



ikQUZ 

I IT C{^y {2/'K)konZ cos Omax) " a/ {2/Tl)konZ COS 

"'''2 (fconz)3/2 



iS{^J (2/Tx)kQnz 



cos Or 



. f^C{^{2/n)konz) - ^J(2/n)kQnz + iS{^{2/n)konz) 

(konz)^^'^ 



(A31) 
(A32) 



where 



C{z) 
S{z) 



cos (^2^^^ 



sin ( — ) c/t. 



(A33) 
(A34) 
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Furthermore, 



(0,^) 



^ikonz cos a cos^/^ « sin « rfo 



^ikonz _ gifconz cos Qmax _ ]^ 

ikonz ikonz 

3 gifcon^ _ I _ ^^^^^ 3 gifco 



cos ^ Clmax 



n.2 cos amax 



1 — ikonz COS a 



{k^nzf 



max 1/2 

COS Clmax 



^J{2/'K)kQnz 



cosa,! 



3 [^l C{^J (2/ tt) fepraz cos 
2V 2 ((2/7r)A;on^)5/2 

3 /7r5(v/(2/7r)A;onzcosamax) - f ((2/7r)A;o'^zcosQ;max)^''^ 



'^2V2 ((2/7r)A;onz)5/2 



2 V 2 



((2/7r)A;onz)5/2 



^7o" (0,-^) = io,zy 



In particular, 



gT\0,0) = -(1-cos 



Hence, along the optical axis, the electric field is 



5/2 



Or 



E(0,0,^) 

and the electric energy density is 



innf 



io,z)+9r (0,2) 



|E(0,0,2)p = 1^.(0, 0,^) 



In the focal point: 



|E(0)| 



2 _ ^'n'P 



^0 



16 2 



■(o,^)+^o" (o,^)r- 



3/2 2 5/2 

TTT 7: COS Q^max T COS Q^max 
15 3 5 



Now 



gl'\0,z)+gl'\0,z) 



Jkonzcosaf 1/2 



cos^^"' a + cos'^^^ a) sin a da 



^ikonz cos a 



sin a da + X(2;) 



(A35) 
(A36) 



(A37) 
(A38) 

(A39) 



(A40) 



(A41) 



(A42) 



44 



where 

/"Omax 

I{z) = / e'*^o""™^"(cos^/2 ^ + pog3/2 ^ _ 2) sin a da. (A43) 
Jo 

For small numerical aperture: 



max 



2 



4 ^max 

(A44) 



where 0{a^^^) denotes a term which modulus is smaller than Ca^^^, where C is a positive 
number that is independent of z. Furthermore, T{z) can be estimated as follows 



/•«max 

< / (2 — cos^^^ a — cos'^^^ a) sin a da 
Jo 



2 2 

2(1 - COSttmax) " "(1 " COS'^^^ ttmax) " "(1 - COS^^^ ttmax) = O(a^ax)- (A45) 



Hence, 

"sin (^aLx)'' 



kpnz 2 
4 max 



,4 I /n/'^,6 



O(aLx)- (A46) 



The first term is identical to the intensity distribution along the optical axis in the scalar 
paraxial theory. This shows that when the numerical aperture is small, the focal depth of 
the intensity in the vectorial theory is the same as in the scalar paraxial theory. 
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